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Abstract 

We study the mean field dynamics of a model introduced by Phan et al [Wehia, 2005] of a 
polymorphic social community. The individuals may choose between three strategies: either not 
to join the community or, in the case of joining it, to cooperate or to behave as a free-rider. Indi- 
viduals' preferences have an idiosyncratic component and a social component. Cooperators bear 
a fixed cost whereas free-riders support a cost proportional to the number of cooperators. We 
study the dynamics of this model analytically in the mean field approximation for both parallel 
and sequential updating. As we vary one of the parameters while keeping the other parameters 
fixed, the phase diagram experiences a rich class of bifurcations. Noticeably, a limit cycle is 
shown to exist in both parallel and sequential updating, under certain parameter settings. A 
comparison of the analytical predictions with computer simulations is also included. 

Keywords: Social Networks Interactions, Dynamical Systems, Bifurcation. 



1 Introduction 



Recently 0H] a model was proposed to analyze social organizations whose members are 
expected to cooperate to a public good. Basic evidence on several kinds of communities, 
as well as data obtained in public goods experiments [31 @] , reveal a rough partition 
between individuals that cooperate to the public good and pure consumers (also called 
free-riders). This polymorphic configuration seems to be a stable form of organization. 

The model corresponds to the following situation: individuals have to decide 
whether to join an organization or community. The members of the organization have 
to contribute to a task whose realization is beneficial for everybody. Cooperators bear 
a fixed cost for producing the public good for the community. The surplus of all the 
individuals, cooperators or not, increases proportionally to the number of cooperators. 
Individuals that do not cooperate are punished by cooperators through costless moral 
disapproval, which may be either a subjective moral burden, or a true sanction. This 
cost, proportional to the fraction of cooperators, is idiosyncratically weighted. 

In the present paper, we consider the dynamical evolution of the model analyti- 
cally, within a mean field approximation. Parallel and sequential updating are studied 
through the evolution of the corresponding map and flow respectively. They present 
different behaviors, depending on the models parameters. 

The most striking result of our analysis is that for some range of parameter values 
and initial conditions, the trajectories of the flow reach the fixed points determined 
by Phan et al. through quasi-cyclic paths that may last for very long times. The 
corresponding map exhibits a limit cycle. Thus, under such conditions, macroscopic 
fractions of individuals change their strategies successively. In experimental economics 
settings that correspond to parallel updating, true cycles might appear. In sequential 
updating, the transient oscillatory behavior of the system may be very long lasting, 
and should thus be observable in actual systems. 

The paper is organized as follows: next section presents the details of the model 
and its equilibrium fixed points. Section |3] is devoted to the study of a mean field 
approximation of the dynamical equations. These describe the evolution of the fraction 
of cooperators and of free-riders. We consider both parallel and sequential updating, 
and study the corresponding map and flow. We show that the system may exhibit 
strong oscillations of the fraction of cooperators and of free-riders, and that limit 
cycles may exist for some range of the parameters. Numerical simulations presented in 
section |3] show that the predicted behaviours exist in finite size systems and might be 
observable in actual systems. The paper ends with a discussion and some conclusions. 



The basic economic model analyzed in this paper was introduced by Phan et al. [Q. It 
considers a system of ./V agents that must choose one among the following strategies: 
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(to join the community and cooperate) 
(to join the community and free-ride) 
(not to join the community) 
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Each individual i has a private idiosyncratic willingness to join the community, hi. 
Its mean value over the population is h and we write: h{ = h + yi where yi is a 
quenched random variable with zero mean. Individuals that join the community have 
a social benefit proportional to the fraction of the population that join the community, 
weighted by a coefficient j, and an additional payoff proportional to the fraction of 
cooperators, weighted by a constant g. Then, depending on whether they cooperate 
or not, they bear different costs. A cooperator bears a fixed cost c. Free-riders do not 
bear this cost, but instead support a moral punishment proportional to the fraction 
of cooperators. This punishment is weighted by an idiosyncratic positive constant x%. 
All the parameters of the model (h,j,g,c as well as the values of {yi} and {x{\ for 
1 < i < N) are measured in units of the variance of the random variable yi. 
The fraction of cooperators is: 

1 N 
i=i 

where 5 denotes the Kroenecker delta. That of free-riders is: 

1 N 

i=l 

so that the fraction of individuals that join the community is rj a = rj c + rjf < 1. 
Each agent chooses the strategy that maximizes his surplus function, 

Vi (Si) = (A, + Bi) 6 Si;1 + Ai 6 Si ,2- (1) 

Ai is the surplus of joining the community being free-rider, and Bi is the bonus of 
cooperating: 

Ai = h + yi + j(r/ c + rjf) + (g - Xi)r] c , (2a) 

Bi = XiT] c - c. (2b) 

From equation^), the best response of agent i to the neighborhood's behavior can be 
rewritten as follows: 

S { = 1, <^> Ai + Bi > and B { > 0; (3a) 
S { = 2, ^> Ai > and B { < 0; (3b) 
Si = 3, otherwise (3c) 

2.1 The Cluster Structure 

In order to get some insight about the problem, it is useful to consider a plane whose 
axes are the quenched random variables x (abscissas) and y (ordinates), as in figure 
n Each individual i is represented by a point according to his values (xi,yi). The 
lines represent the marginal individuals whose utilities are at the boundaries between 
different optimal strategies. These are: 

V = ym = -h-j(r] c + r] f )-grj c + c, (4) 

V = ym-c + xr] c , (5) 



Then, according to (f5|). individuals in region I will choose to cooperate, in region II to 
free-ride and in region III not to join the community. 



~B ^ — . x — x y, 
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Ai + Bi = 0; y = y n 
Ai = ; y = y m - c + xrj m 
III 



Figure 1: Boundaries between regions corresponding to the different individual choices, 
as a function of the values of the quenched random variables. 



2.2 Fraction of cooperators and free-riders at equilibrium 

The fixed points of the system for different distributions of the quenched random 
variables {xi,yi} have been discussed in and [2] in the mean field approximation 
N — > oo. We summarize them here for completeness. 
Introducing the complementary cumulative functions: 

/•oo 

G X (C) = 1-F X (C) = J f x (m (7) 

where x £ {^^}) an d F\(C) is the cumulative distribution of f x , rj c and rjf have the 
following expressions in terms of the G functions and the marginal values Q) : 

Vc = G x (x m ) G Y (y m ), (8a) 

r lf= fx(x) G Y (y m - c + xrjc) dx, . (8b) 
Jo 

The solutions to equations (jHJ) give the fraction of cooperators and free-riders at equi- 
librium as a function of the parameters of the problem, namely h, j, g, c, d (these pa- 
rameters are measured in units of the variance of the distribution of y). Notice that 
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r] c = is always solution of (|8a[) . in which case all the members are free-riders, and the 
problem reduces to that of the simple social interactions model considered in [SJ El an d 
references therein. It is important to stress that, due to the non symmetric interaction 
Xi, there is no reason that the stationary states of the system be exclusively fixed 
points, as is the case for symmetric interactions 1 . 

In this paper we consider the most interesting of the cases analyzed in [I] , with the 
Xi following a uniform distribution of finite width d: 

f x ( x ) = 1 for < x < d, (9a) 
fx(x) = otherwise. (9b) 
and the yi distributed according to the following probability density function: 

- iseW (10) 

The cumulative function corresponding to (|1()|) is the logistic distribution, and its 
complementary function is: 

GV(z) = l/[l + exp(z)]. (11) 

In this case, following Phan et al. pQ, equations © may be written as follows: 

i lc = (l-^)G Y {c-Z) (12a) 



c 



Vf = — log[l + (e c - l)G Y (c - Z)}, (12b) 

CT] C 

where p = c/d and Z = h + j(r] c + rjf) + grj c . 

Calling G the value taken by Gy(c — Z), one has rj c = (1 — ^) G, which gives r\ c 
(if > 0) in terms of G: 

r ]c = vt[G]^\G{l±[l-^} 1 / 2 } (13) 

Later we'll see that both the + and the — branches may give stable equilibria, though 
the parameter range for which a stable equilibrium exists for the — branch is much 
narrower than for the + branch. The function is very close to its asymptote, 

G — p, for all the values of G. From (J13|) one gets that if r\ c ^ 0, then r\ c > 2p. Equality 
can occur when j = g = 0. 

Equation (|12b|) can also be parameterized in terms of G by defining 



_ j(vc + y f ) + gyc (u) 

3 + 9 

Introducing (|T3|) and (|12b|) into (|Ti|) we obtain an equation for rj: 

V = Vi[G] = vtlG] + J- -4?^ ^log[l + (e c - l)G] (15) 
l±9r]c[G] d 



x In Ising systems, with binary microscopic states and symmetric interactions, the only attractors 
are fixed points for sequential dynamics, and either fixed points or cycles of order two for parallel 
dynamics. 
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Inverting Ijlljl for z = c — h — (j + 5)77, we obtain 77 in terms of G: 

V = m[G] = ~r~~ {c - h - log -}. (16) 
3+9 G 

The possible solutions 77 are then obtained by the intersects of the curves 771 [G] and 
r)2[G]. Introducing the corresponding value of G into (|13f) allows to determine 77 c , and 
then 77/ is deduced by introducing the values of rj c and rj into (|14[) . An example of 
curves r?i[G] and 772 [G] (G G [0,1]) is shown on Figure (J2J), which presents up to 3 
intersects (for h — c = —3.3), although one can check that at most 2 correspond to 
stable equilibria. Notice that this analysis allows to determine the fixed points of the 
system, but doesn't give any hint about the existence of cycles. 




Figure 2: Solution by curve intersection - a case of interest. 



3 Mean field dynamics 

We are interested in the temporal evolution of the system within a repeated game 
setting, in which individuals have to choose their best strategies Si(t) based on the 
available information. We assume that each time an agent has to make a decision, he 
has the exact information of the global proportions of cooperators and free-riders at 
the preceding outcome, rj c (t — 1) and r]f(t — l), and uses these quantities to estimate his 
utility ((21). Such a dynamics is called Cournot best reply in economics literature. In 
simulations, starting from initial guesses ?7c(0) and 77^(0), the updating is said to be in 
parallel if all the individuals in the population first determine their best strategies based 
on the preceding outcome, and make their decisions simultaneously afterwards. At the 
opposite, in random sequential updating, a single individual is selected at random and 
asked to make his decision at each time step. The latter dynamics simulates systems 
where the individuals make their decisions without any time coherence. In order to 
compare the time scales of both dynamics, it is usual to consider that iV sequential 
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time steps are equivalent to one parallel update. Intermediate updating schemes may 
be implemented, but we only consider these two extreme cases in this paper. 

Referring back to Figure since the boundary lines depend on the values of rj c (t) 
and rjf(t), they will shift in the course of updating according to the perceived pro- 
portions of cooperators and free-riders. Individuals whose values of yi and Xi lie close 
to the boundaries are susceptible to small changes of rj c and rjf. their strategies may 
change with time, and in turn induce changes in those of the others. 

We are interested in the way the system reaches its stable states upon successive 
updates, that is, the path followed by a point representative of the system's state in 
the plane (r] c ,r]f). Such paths should either end up in one of the stable fixed points 
determined by the previous analysis, or get trapped in other types of attractors if they 
exist. In the following we consider separately the two updating schemes, since the 
corresponding dynamic equations are different. 



3.1 Parallel dynamics: two-dimensional map 

In parallel dynamics, one assumes that the agents are updated simultaneously start- 
ing from an arbitrary initial configuration, i.e. the perceived r\ c and rjf used by all 
the agents to estimate their utilities are the same. The boundary lines between the 
three regions of figure ^ which depend on rj c (t) and rjf(t), partition the population 
at time t according to their estimated utilities: all the individuals whose idiosyncratic 
parameters lie in region I will choose to cooperate, those in region II to free-ride 
and those in region III not to join. In this way, using equations ©, the dynamics 
can be formulated as the following two-dimensional deterministic map in the simplex 
S := {{Vc Vf) ■ Vc > 0, rjf > 0, rjc + Vf < 1}: 

n : iv'c Vf) = (p(Vc Vf),q(Vc rjf)) ( 17a ) 

pOO pOO f'Xm POO 

= ( / fx(x)dx / f Y (y)dy, / / fx(x)f Y (y)dydx). (17b) 

Jxm Jy m J -oo Jym-c+xric 



In the case considered in this paper, the density functions are © and (flflj) . Since fx 
has a bounded support, the first equation vanishes if < rj c (t) < p. Thus, if the initial 
value of rj c is such that r] c (0) < p, then rj c (l) = and a state with no cooperator 
is reached after a single parallel update. Afterwards, on the axis rj c = 0, rjf evolves 
according to 



Vf(t + l) = / fy(y)dy. (18) 

J-h-jr)f(t) 

The equilibrium value of rjf is obtained replacing rjf (t) by rjf in the above equation. 
It has been shown [7] that there is a critical value of j, js, such that for j < js, there 
is a single equilibrium. In that case all the points rj c (0) < p will eventually be mapped 
to it. If j > jB equation ((TH|) with rjf{t) = rjf has three solutions. One of them, rjf u is 
an unstable fixed point separating the basin of attraction of the two others, that are 
stable. To summarize, when r/ c (0) < p, after a first time step the system is mapped to 
the axis rj c = and then evolves, following (|18|) . either to one fixed point (if j < jb) or 
to one of two fixed points (if j > js) depending on whether r]f(0) > rjf u or r//(0) < rjf u . 
In the case of the logistic distribution considered here, jb = 4. 
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When ?7 C (0) > p, the dynamics is most fruitfully studied numerically. In two- 
dimensional maps or flows, a fixed point is called a sink if all the points in its neigh- 
borhood converge to it, a source if these points diverge from it and a saddle if in one 
direction the map or flow converges, while in the other directions it diverges. Together 
these three types are called hyperbolic, which is the only kind of equilibria possibly 
encountered in a structurally stable system, i.e. a system stable with respect to small 
variations of the parameters (in our case h, j, g, c and d). For our two-dimensional 
map, the nature of an equilibrium point (rjcOiVfo) 1S determined by the following Ja- 
cobian matrix 

[ 9q(r]co,Vfo)/dVc (9<?(r/ c o,?7/o)/<% 

Denoting the two eigenvalues of J m as a\ and 0,2, then the fixed point is a sink if \a\ j < 1 
and |«2 1 < 1, a source if |ai| > 1 and | a-2 1 > 1 and a saddle if (|ai| — 1 ) ( | a-2 1 — 1) < 
( [HJ §1-4) . There are two types of curves emanating from each saddle point xq, namely 
the stable manifold, defined as 

W s (x ) = {x G R 2 \Q n (x) -> x as n -> oo} (19) 

and the unstable manifold, defined as 

W u (x ) = {x G M 2 |fT n (x) -> x as n -» oo} (20) 

Q,^ 1 denoting a backward iterate of the map. The boundaries of the basins of attraction 
are usually formed by stable manifolds. Moreover, if a transversal intersection between 
the stable and unstable manifolds of the same saddle point exists, the map will exhibit 
chaotic behavior. 

In view of their crucial importance, we have numerically computed the stable and 
unstable manifolds for each saddle point, with the same set of j, g, c and d as in 
Figure while changing h — c from —4 to —3 in steps of 0.02. The algorithms used 
to compute the stable and unstable manifolds can be found in jH] and Chapter 10 
of ^U|, respectively. The results classified according to the qualitative dynamical 
features are shown in Figures (|3I8[) . Hereafter we provide the increasing values of 
h — c, denoted by Aj, at which the successive bifurcations appear, together with the 
qualitative descriptions of their nature. 

• Ai = —4.0971. For h — c < Ai there exists no intersection between 171(G) and 
172(G) and all the points in S are eventually mapped to the trivial fixed point, 
which is always a solution of the equations. At h — c = Xi, 172(G) begins to 
intersect 171(G) at two points, one of which is a saddle (HJ) and the other is 
a source (X— ) (Figure OJ) ■ The process of simultaneous creation or elimination 
of a pair of equilibrium points is called a saddle-node bifurcation. Note that 
this bifurcation does not alter the overall behavior of the system, since no new 
attractor is created. Thus, only the fixed points rj c = 0, rjf > exist, and since 
j > js, the system may flow to either of the two fixed points for 17 f depending 
on the initial conditions. 
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Figure 3: Map. The stable and unstable manifolds for h — c = —3.7. Throughout this 
article, stable manifolds are represented by dashed lines and unstable manifolds by 
dotted lines. It's well understood from their definitions that an arbitrary initial point 
on the stable manifold will be mapped towards the saddle from which it emanates, 
while one on the unstable manifold will be mapped away from it. Sinks, sources and 
saddles are represented by O, X, and H respectively, and the subscripts describe the 
branch (rj^ (G) /rf[ (G)) on which the equilibrium is found. 
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Figure 5: Map. A separatrix bifurcation between the stable manifold of iJ_ and 
the unstable manifold of H + has introduced a first order transition in the basin of 
attraction. 




P 0.2 0.4 0.6 0.8 1 
'c 



Figure 6: Map. A separatrix bifurcation between the stable and unstable manifolds of 
H + has introduced a stable cycle C s with its own basin of attraction. 
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• A2 = —3.4201. Another saddle-node bifurcation happens on rjf(G) at h — c = A2, 
with a saddle (-H+) and a sink (0+) created simultaneously (Figure Now the 
stable manifolds of H + and H- possess a common end emanating from X + , and 
together they divide S into two regions. Points outside the region bounded by 
the stable manifolds should converge to the trivial fixed point r/ c = 0, r/f > as 
before, and those inside to O+. 

• A3 = —3.3617. For A2 < h — c < A3, the unstable manifold of H + approaches the 
stable manifold of until they coincide at h — c = A3 (the common manifold 
is called a separatrix between H + and H_). Beyond this value, the unstable 
manifold of H + folds back and converges to 0+ while the stable manifold of 
H- hits the simplex boundary r/f = (Figure EJ). Now the basin boundary is 
determined solely by the stable manifold of H^, causing the basin of attraction 
of + to experience a sudden expansion. 

• A4 = —3.3504. For A3 < h — c < A4, the stable and unstable manifolds of H + 
approach each other until they coincide at h — c = A4. Beyond this value, the 
former ends up at the boundary r/f = of S and the latter folds back into a 
stable limit cycle C s around X + (Figure EJ). At this transition C s is introduced 
into the system as a new attractor with its own basin of attraction delimited by 
the stable manifold of H + . 

• A5 = —3.2726. For A4 < h — c < A5, C s shrinks and eventually merges with 
X + at h — c = A5 to produce a sink 0' + (Figure E| • The process of transition 
between a sink and a source with the simultaneous appearance or disappearance 
of a limit cycle is called a Hopf bifurcation. 

• A6 = —3.1146. For A5 < h — c < Xq, 0' + and H + approach each other until 
they are annihilated at h — c = \q. We are left with and + , the stable 
manifold of serving as the basin boundary (Figure EJ). This topology persists 
for arbitrarily larger values of h — c. 

Summarizing the analysis above, we can identify three types of transition: 

1. saddle- node bifurcation; 

2. Hopf bifurcation; 

3. separatrix bifurcation between two saddle points, or one saddle point with itself. 

Type 1 can be directly determined from the number of intersections between r\\ (G) and 
772(G) (cf. figure EJ), or from the fact that J m has an eigenvalue +1. Due to the lack of 
symmetry in our system, transcritical and pitchfork bifurcations never occur ([S] §3.4). 
In case that J m has an eigenvalue —1, one can have a period-doubling bifurcation ([B] 
§3.5). This never happens in our map, since one can easily verify that 

J m (l, 1) > 0, J m (2, 2) > 0, A(J m ) > 0, (21) 

which excludes the possibility that J m has a negative eigenvalue. 
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Types 2 and 3 cannot be observed from the intersections between 771(G) and 772(G) 
and must be determined numerically. Type 2 remains a local bifurcation with J m 
having a pair of conjugate complex eigenvalues with unit modulus. In contrast, Type 
3 is a global bifurcation with the area of the basin of attraction experiencing a sudden 
jump or, in physicist's terms, a first-order transition. Numerically we have found no 
evidence for homoclinic intersection leading to chaos, but the question as to whether 
chaos exists in our two-dimensional map remains open. 

3.2 Sequential dynamics: two-dimensional flow 

In sequential dynamics, one updates one agent chosen at random in each time step. 
It is usual to consider N (the number of agents in the system) successive random 
sequential updates (called one Monte-Carlo step) as being comparable to one step of 
parallel updating. 

In each individual update, the expected displacement of the pair (r] c ,rif) is 

(Arj c ,Arjf) = j^(p(VcVf) ~ Vc,q(Vc,r)f) - Vf), ( 22 ) 

where p and q are defined in (|17|) , Taking the continuous limit N —> 00, we have the 
following set of differential equations 

dri c /dt = p(r) c ,r]f) - r] c , (23a) 
drjf/dt = q(r) c , 77/) - rj f , (23b) 

with the time unit being one Monte-Carlo step. Now the community evolves as an 
autonomous system with planar phase space, or a two-dimensional flow in S. 

The three generic types of equilibria, namely source, sink and saddle exist in a 
structurally stable two-dimensional flow as well. However, the nature of an equilibrium 
point (rjcOiVfo) is now determined by the Jacobian matrix of the flow 

dp{r]cO:Vfo)/dVc - 1 dp(r] c0 ,r] f0 )/dr] f 
09(»7co,»7/o)/&fc ^(^0,77/0)/% - 1 / 

Denoting the eigenvalues of J/ as a\ and 02, then the point is a sink if a\ < and 
C12 < 0, a source if a\ > and > and a saddle if a\ai < ([Sj §1.2~1.3). Note 
that the condition for an equilibrium to be a saddle gives the same inequality as in the 
parallel case. Therefore, the set of saddle points in sequential updating coincides with 
that in parallel updating. However, the set of sinks or sources are not necessarily the 
same in both cases, since the governing inequalities are quite different. The definitions 
of the stable and unstable manifolds are analog to those in the two-dimensional map 
with the discrete time step replaced by the continuous time variable t. However, 
we note that the two-dimensional flow system, according to the Poincare-Bendixson 
Theorem (0 Theorem 1.8.1), will never go into chaos. 

We have computed the stable and unstable manifolds as for the map. The results 
shown in Figures (|9I13|) . Hereafter we provide the bifurcation values of h — c, denoted 
by fn, together with the qualitative descriptions of the nature of the bifurcation: 
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Figure 9: Flow. The stable and unstable manifolds for the flow. Compare with Figure 

El 
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Figure 10: Flow. A sink O together with an unstable limit cycle C u has emerged 
from the source X_ through a Hopf bifurcation. 
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Figure 11: Flow. A separatrix bifurcation has eliminated C u and expanded the basin 
of attraction of the sink. 




Figure 12: Flow. H + and O+ have been created from a saddle node bifurcation on 
r)f(G). Compare with Figure El 



15 



Hi = \\ = —4.0971. As in the parallel case, a saddle node bifurcation occurs at 
h — c = fj,i and produces a saddle H- and a source X_ (Figure 



• \il = —3.6172. A Hopf bifurcation occurs at h — c = fi2 and produces an unstable 
limit cycle C u together with a sink 0_ from the source A_ (Figure I1U|) . The 
basin of attraction of 0_ is the area surrounded by C u . Note that we've been able 
to choose h — c such that the sink remains in the negative branch, which implies 
that variation of rj with respect to h — c does not always qualify to determine the 
stability of the fixed point. 

• ^3 = —3.4980. A separatrix bifurcation of the stable and unstable manifolds of 
H- eliminates C u (Figure ITT|) . 

• m = A2 = —3.4201. A saddle node bifurcation occurs at h— c = \x± and produces 
a sink + and a saddle H + (Figure IT2*|) . 

• 115 = Xq = —3.1146. A saddle node bifurcation occurs at h — c = ^5 and 
annihilates the saddle H + and the source 0' + . The topology persists for larger 
h — c ( Figure IT3|) . 

Thus, we have the same types of bifurcation as in the parallel case, and Peixoto's 
Theorem (jHJ Theorem 1.9.1) indicates that these are the only types of bifurcation ex- 
pected in a two-dimensional flow system. However, let's note that the overall topology 
of the dynamics depends strongly on their order of occurrence. For example, if the 
separatrix bifurcation happens before the Hopf bifurcation, a stable limit cycle can be 
created just like in the parallel case (Figure El- 




P 0.2 0.4 0.6 0.8 1 



Ti 
'c 

Figure 13: Flow. A final saddle node bifurcation has annihilated 0' + and H + . Compare 
with Figure El 
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4 Comparison with Computer Simulations 

We simulated systems of N = 1000 agents, corresponding to the parameter values of 
figures fT5l to I20 ( with parallel dynamics. Starting from different initial conditions, we 
determined the trajectories in the phase space (r] c , rjf). The results, presented on figures 
El to 1201 show an excellent egreement with the analytical results. The complexity of 
the map diagram is reflected on the winding trajectories. One of the most striking 
results of the simulations is that the fraction of cooperators and free-riders oscillate in 
time, even in systems for which do not expect to have cycles. 

5 Conclusion 

In this paper we have studied the time evolution of a social system with cooperators 
and free-riders. We analyzed the two-dimensional dynamical system and obtained 
different phase diagrams corresponding to various parameter settings. The fact that 
a rich class of bifurcations can occur within the narrow tubular region enclosed by 
the curves h — c = — 3.7 and h — c = — 3.1 in Figure (J2J) is in itself truly remarkable. 
Physically this phenomenon has its root in the strong asymmetry inherent in the 
interaction between agents, due to the idiosyncratic weights that the free-riders give 
to the social disapproval. Moreover, in the context of parallel and sequential updating, 
we have sometimes distinctly different phase diagram for the same parameter setting. 
A natural question to ask is what happens in the intermediate case, i.e. some agents' 
decisions are taken simultaneously while others' decisions are taken independently of 
each other. Numerical simulations show that the analytical results obtained in the 
thermodynamic limit are valid for systems of N = 1000 individuals. The dynamical 
behaviour of the simulated systems show winding trajectories in very large regions of 
the phase diagram. Oscillations of the fraction of cooperators and free-riders are thus 
expected in such systems. We are currently performing simulations with sequential 
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Figure 14: Flow. A stable limit cycle in the two-dimensional flow. 
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Figure 15: Examples of trajectories of a simulated system, starting from different 
initial conditions, corresponding to the map of figure El 




Figure 16: Examples of trajectories of a simulated system, starting from different 
initial conditions, corresponding to the map of figure HJ 
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Figure 17: Examples of trajectories of a simulated system, starting from different 
initial conditions, corresponding to the map of figure [5J 
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Figure 18: Examples of trajectories of a simulated system, starting from different 
initial conditions, corresponding to the map of figure 
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Figure 19: Examples of trajectories of a simulated system, starting from different 
initial conditions, corresponding to the map of figured 
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dynamics. Already obtained results also exhibit the above mentionned oscillatory 
behaviour. Future work can also focus on the dynamics in the case of non-global 
neighborhood, which probably requires more sophisticated techniques from the theory 
of dynamical systems. 
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Figure 20: Examples of trajectories of a simulated system, starting from different 
initial conditions, corresponding to the map of figure [Sj 
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